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Abstract 

The technique of vector differentiation is applied to the problem of the derivation of multipole 
expansions in four-dimensional space. Explicit expressions for the multipole expansion of the 
function r n Cj(r) with r = ri + T2 are given in terms of tensor products of two hyperspherical 
harmonics depending on the unit vectors fi and ?2. The multipole decomposition of the function 
(ri • T2) n is also derived. The proposed method can be easily generalised to the case of the space 
with dimensionality larger than four. Several explicit expressions for the four-dimensional Clebsch- 
Gordan coefficients with particular values of parameters are presented in the closed form. 
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I. INTRODUCTION 



The quantum angular momentum theory has proved to be a very efficient tool in many 
physical problems. Among the basis objects of that theory are the spherical harmonics which 
are the solution of the angular part of the Laplace equation. For spaces with dimensionality 
larger than three, the eigenf unctions of the angular part of the Laplacian are often called 
"hyperspherical harmonics" (HSH). From the point of view of quantum mechanics HSH are 
eigenfunctions of the total angular momentum operator. That explains why HSH have found 
wide use in the theory of quantum few-particle systems. 

The theory of three-dimensional spherical harmonics is well developed since 
J.C Maxwell's work on electromagnetic theory. In particular, the use of spherical har- 
monics allows one to simplify the calculation of various angular integrals. However, before 
such an integral can be calculated, one has to re-write the integrand functions in terms of 
spherical harmonics. The decomposition of a some function in terms of spherical harmonics 
is called the multipole expansion. There is a large number of known multipole expansion 
formulas for different functions in three-dimensional space, see e.g. 

The expansion formulas in many-dimensional space are less known. The multipole ex- 
pansion for the exponential function exp(ia • b) has been derived in [3j]. This exponential 
function is a scalar and its multipole decomposition is equivalent to the expansion in terms 
of Gegenbauer polynomials. In the paper |3| also the differential method for deriving the 
multipole expansions of scalar functions of the kind /(|ri — r2|) has been developed. 

Generally, the explicit expressions for the coefficients of multipole expansions (so-called 
"multipole coefficients") can be obtained using the orthogonality properties of spherical 
harmonics. This means that the multipole coefficients are determined by calculating the 
overlap integrals between the function to be expanded and the corresponding HSH. It is 
clear that in spaces with higher dimensions those overlap integrals will become more and 
more complicated. 

The main goal of the present paper is to develop the differential technique for the deriva- 
tion of multipole expansions. In the three-dimensional space such a technique is already 
known (^J . In the present paper, the multipole expansions in four- dimensional space are con- 
sidered. The study of multipole expansions in four-dimensional space is motivated mainly 
by the fact that four-dimensional spherical harmonics represent the wave function of the 
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hydrogen atom in momentum space M. Such harmonics are also used as Sturmian basis 

HQ 

functions in many-centre Coulomb problems, see e.g. jal3|- We note also that the classical 
electrodynamics is built in the four- dimensional Minkowski space. Formally, this space cor- 
responds to the space of vectors whose three (space) components are purely imaginary and 
the fourth (time) component is real. The rotations in Minkowski space are represented by 
the Lorentz group. The multipole expansion technique developed below remains the same 
also in Minkowski space. 

As an example of the developed formalism, several expansion formulas are derived for 
tensor functions depending on the linear vector combination ri + r2. The proposed method 
can easily be generalised to the case of many- dimensional spaces. However, this requires the 
extensive study of the properties of corresponding Clebsch-Gordan coefficients which is not 
the subject of the present paper. 

The paper is organised as follows. In Section |H] the general properties of four- dimensional 
HSH are discussed. Unlike three-dimensional spherical harmonics where there is only one 
set of their indices (and arguments) exist, HSH in four-dimensions can be parametrised in 
several ways. These are analysed in ref. Pj. Below we consider only two sets most important 
for practical applications. Keeping in mind the hydrogenic 0(4)-symmetry, we denote them 
as spherical-type and parabolic- type HSH. In Section HTT1 the multipole expansion differential 
technique for four-dimensions HSH is developed. It is based on the rearrangement formula 
for the exponent exp(a- V) where V is the gradient operator and a is an arbitrary vector in 
four- dimensional space. Although the expressions for the Clebsch-Gordan coefficients (CGC) 
of the 0(4) group are known (see e.g. |illlf|). the definitions used in common literature may 
differ in phase. Therefore, below the explicit expressions for CGC are considered. Analysis 
shows that for some particular values of indices, CGC can be written in a compact closed 
form. The corresponding expressions are given in Appendix [X] 



II. HYPERSPHERICAL HARMONICS IN FOUR-DIMENSIONAL SPACE 

In this section the properties of irreducible tensors in four-dimensional space are analysed. 
The Cartesian components of the vector r in four-dimensional space we denote as (x, y, z, zq). 
We will also use the set of hyperspherical coordinates r, 8q, 8, <p which are connected to 
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Cartesian coordinates by means of identities 

x =r sin 9 sin 9 cos (j), 
y =r sin 9n sin 9 sin d>, 

(1) 

z =r sin 9 cos 9, 
z =r cos# . 

As is seen 9q is the angle between the ^o-axis of the coordinate frame and the vector r. 
(We note that in Minkowski space the length is defined as r 2 = z 2 — x 2 — y 2 — z 2 and the 
replacements 9 — > i9 should be made.) The inverse relations are 

r 2 = x 2 + y 2 + z 2 + z 2 , 

9 = arccos — , 
r 



y/x 2 + y 
9 = arctan — 



(2) 



4> = arctan — . 

x 

(The inverse relations in Minkowski space can be derived similarly.) 

By definition, an irreducible tensor is a traceless^l tensor, symmetric with respect to 
the interchange of any pair of its indices The count of different components of the 

rank-j irreducible tensor in four- dimensional space gives the number (j + l) 2 . This means 
that the irreducible tensor can be labelled with two indices each of which runs from to j. 
Tensor products Cj^ir) will be denoted as parabolic-type set of HSH. There exists another 
possibility of choosing the set of tensor indices. Namely, one can use indices A and a so that 
A varies from to j and a runs from —A to A. The total number of all possible combinations 
of indices A and a remains, of course, equal to (j + l) 2 . The set of functions Cj^ a (t) will 
be referred to as spherical-type HSH. 

It is easy to see that the irreducible tensor product of n vectors r satisfies the Laplace 
equation, 

A{r}, iM)l/ = 0, 

cP_ ( 3 ) 

dx 2 dy 2 dz 2 dz 2 

Indeed, the tensor {r}j^ u is a homogeneous polynomial of the order j with respect to 
the components r\^ v . The action of the scalar Laplace operator A on such a polynomial 



decreases its order by two. However, the result of this action should still be an irreducible 
tensor of the same rank j because the tensor equality must contain only tensors of equal 
ranks. Thus, the above statement is proved by contradiction. We note also that the tensor 
product {r}j ^ u = r J {f }j )fi>v satisfies the boundary condition 

^{r}^ -> 0, (4) 

i.e. it vanishes at the origin. Another tensor solution of the four-dimensional Laplace 
equation which is divergent at the origin has the form r~^~ 2 {r }j,/j,,u- Clearly, the tensor 
product of unit vectors f , {r}j iAt)1/ is independent of the hyperradius r, and, hence, it is the 
function of hyperangles. Therefore, the tensor product {r}j^ u coincides with HSH up to 
normalisation factor. (The same is true also for the set of indices j, A, a.) 

Below, we consider two possible parametrisations for the four-dimensional HSH. In 
sec. Ill Al the parabolic- type set of HSH is introduced. This set is particularly convenient for 
the calculations since the corresponding Clebsch-Gordan coefficients have simple form. The 
spherical-type HSH are discussed in sec. Ill Bl Such harmonics describe the wave functions 
of the hydrogen atom in momentum space labelled with spherical quantum numbers n,l,m. 
The Clebsch-Gordan coefficients for spherical-type HSH have more complicated structure 
comparing to that of the parabolic-type HSH. 



A. Parabolic-type spherical harmonics 

Now let us introduce the hyperspherical components of the vector r according to the 
relations 



From these equations it follows that 



r 2 



(5) 



The square length of r in terms of these components has the form 

i 

E = E(- 1 ) M "" r ** r -*-»' = x2 + y2 + z2 + z °- ^ 



The hyperspherical components of r are connected to the hyperangles as 

r 

r ± i ± i =^p(cos6'o =F isin^ocos^), 



r 

±tM ^\ (8) 
r , i x i = ± i— = sin #0 sin 9e Ttft> 

±2^2 ^ 
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These equation can be re-written in terms of three-dimensional finite rotation matrices U$ v 
parametrised by the direction of rotation and the rotational angle Q], 

r^ = ^=U^(26 ,e, ( f>). (9) 

We define the parabolic-type (or if-harmonics) spherical harmonics Hj tll)V (r) according 
to the relation 

= ^TliW?) = ^fiU${20 o , 9, 0), (10) 

where JjjJ^ are elements of the three-dimensional finite rotation matrix. Written in terms 
of Euler angles these matrix elements are called Wigner's D-functions. 

Let us consider if-harmonics for some specific values of its arguments and indices. The 
evaluation of the harmonics ifj )At)I/ (e z ) depending on the unit vector eo with components 
e = (0,0,0,1) can easily be performed noting that its spherical coordinates are r = 1, 
6*o = 9 = = 0. For this case we have [l|, 

^(6^ = 1^(0,0,0) =(5^. (11) 

Note that eo is directed along the ;?o-axis of the coordinate frame. 
The scalar product of two if-harmonics is defined as 

J/2 

(^■(a) • H 3 (h)) = (-l)^^(a)^-(b) = Cj(a ■ b), (12) 

where Cj-(a-b) is the Gegenbauer polynomial. The expansion formula for the product of two 
hyperspherical harmonics Hj itljV can be derived based on properties of the 3<i-finite rotation 
matrices llj/u, 

31+32 

i=b'i-j2 

where Hj^ iUi .j 2 ^ 2U2 are the "parabolic-type" CGC for the 0(4) group. From the definition 
(110)1 of if-harmonics as 3d finite rotation matrices it follows that H-type CGC are simply 



products of two CGC of the 0(3) group 

rL jl^lv 1 ;j 2 ^2V2 °(jl/ 2 )/tl(j's/ 2 )« (n/2h()V 2 )«' ^ ' 

This expression for CGC directly follows from the definition of spherical harmonics as the 
elements of the finite rotation matrix in 3<i-space. We will not discuss further properties of 
coefficients ff^™,,^,, such as orthogonality etc, since they can be easily deduced from the 
properties of conventional CGC in 3<i-space. 

From the expression (fH^ for CGC one can deduce the triangle rule for the momenta 
jxij2ij- Namely, the momentum j can take the values \ji — j 2 \, \j% — j 2 \ + 2, . . . ,ji + j 2 . 
Thus, the sum j\ + j 2 + J3 is always an even number. In particular, for ji = j 2 we have that 
j = 0, 2, 4 ... , (2ji). 

The tensor product of two spherical harmonics we will denote as {//"^(a) <g> Hj 2 (b)}j and 
it is defined by the equation 

ii/2 J2/2 

{H n {k)®H n {h)} 3 ^ u = Yl E H nT^n^ H n^Amn^^)- (15) 

l*l,"l=—h/ 2 M2,^2 = -i2/2 

Below, we will label the tensor products and spherical harmonics only with their ranks 
unless one needs to derive explicit expressions. The bipolar harmonics of the zero rank are 
connected with the scalar product of spherical harmonics as 

{i?j(a) ® #z(b)} = j^jW*) • Hi(b)). (16) 
For two equal vectors a = b this equation becomes 

{#,(&) ® #,(&)}„ = T^TO*) • fl-j(a)) = I- (17) 

B. Spherical-type spherical harmonics 

The explicit expression for the spherical-type HSH can be derived using the parametri- 
sation of the 3c? finite rotation matrix in terms of the rotation angle and the rotation axis, 

21 2A -I- 1 

ul v (29 , e, 0) = J2H) x ^~y cl; XaX {(^o)c Xa (e, 0), (is) 

A=0 
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where C Aa are modified spherical harmonics connected to the usual 3<i-spherical harmonics 
by C\ a = a/47t/ (2A + l)Y\ a . The generalised characters of the 0(3) rotation group x l \ can 
be expressed in terms of Gegenbauer polynomials []], 



X l x (26 ) = (2A)!! V2TTT J {2 f + A ^^ (sin^^+^cosgo). (19) 

We introduce the set of spherical-type (C-harmonics) HSH C, jAiM by connecting them to 
if-harmonics as 

W f) = V^ £ ^/txA,Av), (20) 

where A = 0, 1, . . . j, and a = —A, —A + 1, ... A. The inverse relation can be derived using 
the orthogonality of 3c? CGC 



A=0 V J 

The explicit expression for C-harmonics can be obtained from (JTSj) and it is 



<W*) = H) x ^l^^ x { /2 (2e )c Xa (0A). (22) 

From this equation follows the formula for complex conjugated harmonics 

C* Xa (*) = (-l) A+a Q, A ,_ a (f). (23) 

The connection of C-harmonics with tensor products of vectors is the same as for H- 
harmonics which is given by eq. (jl(J|) . The scalar product of two C-harmonics is defined 
similarly to the above eq. (fT2|) . 

(C,(a) • C,(b)) = ^q A)a (a)C iAQ (b) = ^(-l) A + a C i|A)a (a)C iiA) _ a (b) = C](a- b). (24) 
We note also the expression for C,(e ) which can be derived from eq. (j2*U|l . 



O ijA , a (e ) = (-0 A y^^xf(0)C Aa (0,0) = VJTT^o^o, (25) 
where we have used the property of the generalised characters |j] 

xf (o) = (j + 1) <V 
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Now we have to establish the expression for CGC for the set of C-harmonics. We define 
those coefficients by considering the expansion of the product of two C-harmonics 

CjiM, ai (r)C j2M , a2 (r) = } ai , ,-,A 9Ct ,Cj,A,a(f ) ■ (26) 

jXa 

CGC in this equation can be calculated by expressing each C-harmonic in terms of H- 
harmonics according to the decomposition (|2T)|) . At this stage the CGC for if-harmonics 
will occur. However, these are already known, see eq. ()14|) . After somewhat lengthy manip- 
ulations with 3d CGC one arrives at the expression 



f 

21 21 i 

2 2 2 



21 21 2 

2 2 2 



} , (27) 



C&^o. = (J + l)v / (2Ai + l)(2A 2 + l)C A ^ lA2Q2 v 

Ai A2 A 

v 

where the table in curly brackets is the 9j-coefficient of 3o?-space. The 9j-coefficient with 
two identical rows is non-zero only if the sum of its indices is an even number Note, 
that the sum j\ + j'2 + j is always an even number. Hence, the C-type CGC is non-zero 
only if the combination Ai + A2 + A is an even number. Thus, the parameter A can be 
|Ai — A 2 |, |Ai — A2I + 2, . . . , Ai + A2. 

The bipolar harmonics of C-type are defined similarly to the if -case (jl5|) , 

®C i2 (b)} iiA)Q = C 'X^n^ C nM«A*)C nM ^)- (28) 

Ai,A2,ai,a2 

The expressions f!16|) . (jl7j) for the scalar product of if-harmonics remain valid after the 
replacement H +-> C . 

We have to analyse the properties of the C-type CGC. Clearly, they follow from properties 
of the 3d CGC and 9j-symbols. For example, the orthogonality relation of Ad CGC has the 
form 

^ ( ^iiA 1 a 1 ;i 2 A2a2 ( ^iiAiai;j2A2«2 = $j,j'$X,X'5 a , a > • (29) 
Ai,A 2 ,oi,02 

The exchange symmetry of Ad CGC is 

f3 Xa _ ( 1)31+32+3 (j3 Xa (on) 

jiXiai; jzX2U2 \ > 32X202: jiMoii 

We note that the above identities hold for both, C- and if -types of CGC. We present also 
one more symmetry relation valid for C-type CGC 

/ntf'Aai _ / i \ji+j2+j+X2+ct2 j ^ /nrJiAiai /oi\ 

°jfiAiQi; 32X20:2 v ' Jl + 1 J ^ a;j2 ^ 2-a2 ' 



At some particular values of its parameters the C-type CGC may be written in closed form. 
For example, this happens at ji = where 



(~lj\oi r r r 



For zero projections, C-type CGC also evaluates to a closed form, 



(32) 



r joo _ I j + 1 

Wx00;i 2 00 - \l {ji + 1} {j2 + 1 y 



(33) 



More explicit expressions for C-type CGC are given in Appendix 1X1 

We have to consider also the re-coupling coefficients in 4<i-space. These re-coupling 
coefficients can depend only upon the ranks of tensors but not on their projection indices. 
Therefore, all tensor identities are covariant, i.e. they are equally valid for H- and C-type 
tensor components. As an example, we consider the re-coupling relation 



{{P a ® Qb}c ® {Rd ® S e } f } k = J2( c + !)(/+ + l)(h+ 1) 

gh 



x {{p a ® i?4 9 gi {g b ® 5 e } h } s 



a b c 
d e f 
g h k 



(34) 



/ \ 2 

a fe c 
2 2 2 



del 
2 2 2 

9 h k 
2 2 2 



> 



(35) 



where the table in square brackets is the 9j-symbol in 4<i-space. Its calculation is more 
simple when using if-type CGC. Omitting the details of derivations we present only the 
result, 

a b c 
d e f 
g h k 

As is seen, the irreducible tensor in 4<i-space may be considered as "double-tensors" in 
3<i-space. 

Finally, we note that HSH are orthogonal, 

Cj,xAt)C},x,At)dn= f\m 2 e de [\mede [ v d<f>c jiXta {e ,e,<f>)c; ltXW {e ,e,<i>) 

Jo Jo Jo 



2tt 
3 + 



-5j t f5\ t \iSa,a>, (36) 
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where the integration is performed over the four-dimensional hypersphere S. Often, HSH 
normalised to unity are used 



= 2 ° hxAi) - (37) 



The orthogonality identity holds also for if -type spherical harmonics, 



2tt 



2 
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III. THE MULTIPOLE EXPANSIONS IN HYPERSPHERICAL HARMONICS 

In this section the differential formalism for the derivation of multipole expansions in 
four- dimensional space is developed. The general formulas are derived in Sec. IIII Al In 
Sec. IIII Bl we consider the multipole expansions of the function |a+ r| n C 3 -(b) where b is the 
unit vector directed along the vector sum (a + r) . 

A. The differential multipole expansion formula 

We start with the conventional Taylor expansion formula in many- dimensional space: 

/(a + r) = ^-(a-V)'/(r). (39) 

1=0 

It can be re-written in the simple symbolic form, 

/(a + r) = e( a - v V(r). (40) 

We have to present this formula in a way suitable for the derivation of multipole expansions 
of the function /(r). This can be done similarly to the three-dimensional case considered in 
Pi . Namely, one has to use the multipole expansion for the operator exponent exp (a ■ V). 
Before doing so, we write the multipole expansion of the exponential scalar product in 
four-dimensional space jjj, 

00 -i+i 

e a- r = ^2 l —2(l + l)J l+1 (-iar)Ci(a • f). (41) 



ar 

1=0 



The Bessel function Ji + \{iar) can be written in terms of a hypergeometric function 



3 



as 



[-ix) l+1 



■W-*»0 = 2 '+i(l + l)! ° Fl(/ + 2;3: /4) " (42) 
11 



This formula allows one to re-write the decomposition (|41|) in the form 

00 i i 00 / 

^ = E W o^i( Z + 2 - « 2 ^/4)C/(a • f) = ^ ° F!(Z + 2; aV/4)(Q(a) . {r} ,). (43) 
;=o ' 1=0 

This equation can be used for the derivation of the multipole decomposition for the scalar 
products (a-r) n . Namely, one has to expand the exponent exp (a • r) and the hypergeometric 
function Fi(l + 2; a 2 r 2 /4) into the power series and compare the coefficients at equal powers 
of ar. This leads to the formula 

[n/2] , 

t \n \( ar \ n sr^ n-2k + l ! _ 

fc=0 v 



n-2feV 

2(^ + 1; 



= "■<->* £ (n-OI(n + + 2)ll ' C 'W- C 'W»- (44) 

For small values of n this formula has been checked explicitly. 

Replacing in eq. (J43|) r with V one arrives at the following equation for the Taylor 
expansion (|4U|). 

00 1 

/( a + r ) = E 4m\ o^iC + 2 ; « 2 A/4)(Q(a) • {V},)/(r). (45) 
1=0 

This is the key equation for the derivation of multipole expansions. 

For further consideration it is necessary to calculate the action of the Laplace operator 
on HSH, 

= (?|-*| + ¥) = £<m = -^W). (46) 

Hereafter, we will label HSH only with their ranks and, for the sake of shortness, we will 
omit all projection indices. 

Using the above formula one can prove the identity 

A/(r)Q(f) = Q(r) ( - f(r), (47) 



where f(r) is an arbitrary function depending on r = |r|. This equation can be presented 
in two different compact forms 

A/(r)Q(f) = C^r) ^/ l+ ^f(r) = ^{iV^-^^ f{r). (48) 
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In the important particular case of f(r) = r n these identities lead to the formula 

AVGfr ) - Cm r"' 2 * {n ~ m V + n + 2 y- ] UQ) 
A r U{v) - d{r) r ( n _ t _ 2k )U( l +n + 2 -2k)\V [ ' 

One can re-write this equation in a form free of factorials 

A fc r"Q(f) = Q(f)r- 2fc 2 2fc ( ~ 2 ~ l - H ) , (50) 



\ " / k \ ^ / h 

where (a)& = T(a + k)/T(a) is the Pochhammer symbol. As is seen, this equation is also 
valid in the case of non-integer values of n. 

B. Multipole expansions of translated hyperspherical harmonics 

Below we calculate the multipole expansion of the function / (a+r) = |a+r| n Cj(b), where 
b = (a + r)/|a + r|. According to eq. ([45)1 . in order to calculate the multipole expansion of 
this function, one has to calculate the operator action 

(C,(&) • {V},)/(r) = (G(a) • {V}0r"Q(f). 

From the general symmetry arguments it is clear that this action can be presented as a 
combination of bipolar harmonics 

(G(a) • {V},)r n C,-(f) = r n - l J2 A ii'{ C i(*) ® c i'(*)h> (51) 

v 

where the coefficients Aw are numbers which depend apart of I, V also on j and n but are 
independent on vectors a and r. 

The independence of the coefficients Au> on vectors a and r allows one to simplify their 
calculation. We will calculate Au> for vectors a having zero length: (a • a) = 0. Noting this 
fact and the expression (}2"4"|) for the scalar product of C-harmonics we arrive at the identity 

(G(a)-{Vh) = 2'/ 2 (a-V)'. (52) 

Now one has to evaluate the action of the operator in the rhs of this equation on the product 
r n Cj(r). This can be done using the chain differentiation rule 

(a ■ V)V^(r) = 2^ 2 (a ■ V)V^{r}, = 2^ 2 £ Q [(a • V) V^] (a ■ V)'- fc {r}„ (53) 

fc=0 ^ ' 
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where (,) = l\/(k\(l — k)l) is the binomial coefficient and V-operators in square brackets 
do not act on the outer terms. The term (a • V)^ fe {r}j can be evaluated using the vector 
differentiation technique described in [3| so that 



(a-V)'- fc W 



J! 



-{{*}i-k ® {r^-j+fc}^ 



(54) 



" (j-l + k)\ 

The calculation of operator action in square brackets in ()53|) simplifies because (a • a) = 0. 
After some simple analysis one obtains 

(a ■ V) k r n ' ] = (n - j)(n - j - 2) • • ■ (n - j - 2k + 2) (a • r )V^- 2fc 



{-2f 



(a • r) r 



k r n-j-2k 



(55) 



Now we have to re-write the combination (a • r) fc {{a}/_£. ® { r }j-i+k}j i n terms of HSH. 
This can be achieved using the tensor re-coupling rules. To illustrate this we note that the 
above construction can be presented in the form 

(a ■ if{{a},_ fc ® {r}^ +fc }, = r^ l+2k ({ a } k ■ {?},) {{a},_* ® r};-z+fc}r (56) 

The re-coupling of the tensor products in rhs of this equation yields 

({a} fc ■ {?}*) {{a},_ fc ® = 2^ 2 - k J> + 1)0' + W + !)(/' + 1) 



x 



k k 
Z — k j — I + k j 
I I' j 



{{a} z ®CHr)}„ (57) 



where we have used the auxiliary identity 



{{a} fc <g> {a}i_ fe }g = {a}/5 ig , 

which follows from the fact that (a • a) = 0. The summation index I' in eq. (|57j) takes the 
values of \j — 1\, \ j — l\ + 2, . . . ,j — I + 2k. Thus, the combination I + 1' + j is always an even 
number. 

The four-dimensional re-coupling coefficient in eq. (J57)) is connected with three- 
dimensional 9j-symbols by means of eq. (J35|) . The resulting 9j-coefficient can be evaluated 



14 



in closed form |]J, so that the four- dimensional re-coupling coefficient becomes 



k k 

I — k j — I + k j 
I V j 



-1 



k\ (j-l + k)\ 



(/ + l)!(j + l)!(fc + l)0' + l) 

r (i+±±?i + 2) r + f 



X 



(58) 



r {idoL + fc + 1) r + fc + 2 ) 

The next step is to insert this equation into eq. (J57)l and substitute the result together with 
eq. (f55]) into eq. This leads to the identity 



(a ■ V)'r»C,(f) = r^2'/ 2 £{{a}< ® C,(f )},-(-l ) W ^ 

11 J ' 



_(-l) fe - // 

2^ r/ - 



r (i+M. + 2 ) r + f 



, (59) 



_ (Z-A:)! V 2 /,r(if + Hi)r(if + H2) 

Here, the summation over /c can be performed analytically. Noting also the identity (J52jl the 
above equation evaluates to 

I' + l 

I' 

1 — 9 — i — n. \ i n — n \ 

(60) 



(C,(fi) ■ {V},)r"C,-(f) = r^2</ 2 £(-l)<{Q(a) ® C,,(*)} 



x 



2 — j — 72 ^ f J ~ n 

(j+l-l')/2 



(l+V-j)/2 

where we have used the fact that (— iy +l+l ' = 1. It is important to note that the derived 
equation is valid for arbitrary vectors a (i.e. not only for zero- length vectors). 

According to eq. (J2SJ), we have to act on eq. (plj) with the operator qFi(1 + 2;a 2 A/4). 
Thus, one has to calculate the operator construction 

,2fe 



F 1 (/ + 2;a 2 A/4)r"-^(f) = ^ 



\ k r n - l C v {r). 



k=0 



2 2k k\{l + 2) 



(61) 



Equation (joTIj) allows one to calculate the action of Laplace operators A fc on the product 
r n ~ l Ci'(r). The resulting series leads to the Gauss hypergeometric function, 



r n-l-2k a 2k f_ 2 + l-l'- n 



k=0 



k\(l + 2) k 



l + l'-n 



-2 + 1-1' -n l + l'-n 



1 + 2;-). (62) 
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For the sake of simpler presentation it is convenient to replace the notations by a — > ri and 
r — > r 2 . Now we can write the explicit expression for the multipole expansion of the product 
r n Cj(r), where r = r 1 + r 2 , 

oo 

r n C 3 (v) = E B if(^) {C,(*i) ® C7,'(f a )}i, 
M'=o 

(ri ' r2) " r2 l rj V.( j + 1){ 2 J (i+ ^ )/2 l 2 J (w ,_ i)/2 
f-2+l-V-n l + l'-n , r?\ 

where the summations over Z, /' are performed over all values at which j + 1 — I' — 0, 2, 4, . . . 
and Z + /' > j. Thus, there is only one infinite summation in the above formula. 

Since the combination (/ + /'— j)/2 is a positive integer number, at negative integer values 
of (j —n)/2 the above multipole series are, in fact, finite sums. Indeed, for (n — j) = 2, 4, . . ., 
the second Pochhammer symbols in the expression for B^' vanishes at all values of /, V 
except those with (/ + /' — j) < \j — n|. 

We note also that for r% = i, r 2 = 1 and j = 0, n = 2N, where N can be arbitrary 
integer, eq. (jHHjl reduces to 

(2i (h ■ r 2 )f = E Sf^^, l){Ci(h) ® a(f 2 )}o- 

l=N,N-2,... 

Thus, it is the multipole expansion for the powers of a scalar product which has been derived 
above in a different way, see eq. (jUJ). 

For finite multipole decompositions the question of convergence does not occur. The 
infinite multipole series is convergent only when T\ < r 2 . In the opposite case (i.e. at 
7*1 > t* 2 ) the replacement rj <-> r 2 recovers the convergence of the corresponding series. 

Below we present several explicit examples of the multipole expansion (Jfillj) . There are 
two situations when they have particularly simple form. Namely, at n = j or n = — j — 2 
the product r n Cj(r) satisfies the Laplace equation 

Ar^(f) = A-L<7,(f) = 0. 

As a consequence, the Gauss hypergeometric function in is equal to unity. This is 
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clearly seen also from eq. (|61|). The corresponding multipole expansions are 

r j Ctf) = E (^W^C^i) ® C^^)};, (64) 

= B" 1 )'^ (' + | + X ) {^(fi) ® ^(r 2 )}„ (65) 

where r = i"i + r 2 . 

Let us consider one more example of the multipole expansion eq. (jHSJ)- At j = we have 
to expand the scalar function |ri + r 2 |™. In this case 1 = 1', and noting the identity 

{C,(*i) ® a(f 2 )}o = J^(Ci(r i) • Cj(? 2 )) (66) 
we can write the multipole decomposition of the function \r\ + r 2 | n as 

x {dih) ■ Q(r 2 )). (67) 

We note that for n being an even integer number, the Gauss hypergeometric function reduces 
to the associated Legendre polynomial. As is seen, the derived multipole series are convergent 
only for r\ < r 2 . In the opposite case r\ > r 2 , the replacement r x «-> r 2 should be made in 
the above formula. 



IV. CONCLUSION 

The main results of this paper are the multipole expansions given by eqs. (J44j) . (J63|) . 
dMl), f|63j) and (jBTf) . The derivations given in the paper can be directly generalised on the 
case of space with dimensionality larger than four. According to eq. (JoTj) . this requires the 
knowledge of corresponding re-coupling coefficients. Since the four-dimensional HSH are 
proportional to the wave functions of the Hydrogen atom, various integrals involving those 
functions can be calculated using the technique of multipole expansions developed above. 
Examples of such calculations will be reported in a forthcoming publication. 

Finally, we note that eq. (JBlIj) can be used for the derivation of the general multipole 
expansion of the function f(r)Cj(r) where r = ri + r 2 . Namely, we write the Taylor series 
for the function f(r), 

oo 

f(r)C J (r) = J2fnr n C J (r). (68) 

ra=0 
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Each term of this equation can be decomposed using eq. (JUHj) . which leads to the equation 



/(r)Q(f) = Y, c u\n,r 2 ) {Ci(fO ® C,(f 2 )L-, 



(69) 



where the coefficients Cjjy are defined by 



cg ) (ri,r i ) = X;/^ ) (n,r 2 ). 



n=0 



The expressions for the parameters B^\ri,r 2 ) are given by eq. 



(70) 
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APPENDIX A: EXPRESSIONS FOR FOUR-DIMENSIONAL CLEBSCH- 
GORDAN COEFFICIENTS 

The expression for C-type CGC simplifies for j = ji + j 2 . In this case the 9j-symbol 
in the definition (|27|) of C-type CGC can be re-written in terms of "id CGC, see eq. (9) of 

n 

Sec. 10.8.3 of |l[. As a result, one arrives at the equation 



C 



C?i+J2)Aq 



c 



AO 



c 



Xa 



Jl!j2! 



jiAiai; ji\ia.i w AiOA20 w Aiai Xioli j\\ 

(3i+32 + A + 1)! ( 3l +j 2 - A)! (2A a + 1) (2A 2 + 1) 



1/2 



• (Al) 



LO'i + Ai + 1)! (3i - Ai)! (32 + A 2 + 1)! (j 2 - A 2 )! (2A + 1) J 

This equation simplifies significantly for Ai = ai = 0. In this case the three-dimensional 
CGC are equal to unity, so that 



^0'i+j2)Aa 

°iioo ; j 2 ; 



32\ 



2 A 2 a 2 - °A,A 2 0a,a 2 ^ + 

For A = a = this equation reduces to 



(ji+j2 + A + l)!(ji+j 2 -A)! 
0'i + l)(j2 + A + l)!0' 2 -A)! 



1/2 



(A2) 



r »(jfi+j2)00 
°ji00; j 2 00 



jl + J2 + 1 
(Jl + l) (J2 + 1). 



1/2 



(A3) 
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which is in agreement with the general equation ()33|) . We consider also the particular case 
of eq. (jAl|) for maximal values of projections A, 



l(jl+j2)oi 



(A4) 



For maximal possible values of both A- and a-projections one has 



C 



{31+32) (31+32) {31+32) 



31313^323232 



(A5) 



We present also the expression for CGC with j = j 2 — j\, which can be derived from (|A1|) 
and the symmetry relations (jHOJ) and (|3*T|). 



/^>0'2-il)AQ 



^AOAiO'-'Aiai A 2 a 2 



(J2 + I)! 

0'2 + A 2 + l)!(j 2 -A 2 )!(2Ai + l) 



1/2 



(A6) 



(ji + Ai + 1)! (ji - Ax)! (j 2 - ji + A + 1)! (j 2 - ji - A)!. 
It is also of interest to consider CGC with some set of projections equal to zero, e.g. 
Ai = «i = 0. In this case, A = A 2 and a = a 2 and the 9j-coefficient in the definition (J27)l of 
CGC reduces to the 6j-coefficient jl|, 



^jiOOyzAa 



(_l)*+(h+32+j)/2 j + 1 

1 ' v / 7H~I 



\ 3 3 
A 2 2 

h h n 
222 



(A7) 



In particular, for j\ = 1 we have that j = j% ± 1 and evaluating the 6j-symbol to its explicit 
form we obtain 



r (j-l)\a _ V (j - A) (j + A + 1) r (j+l)\a 
°100;jAa ~~ f • , 1\ /n ' W0;j\a 



(j + l)V2 



y/(j-A + l)(j + A + 2) 
(j + 1)V2 



(A8) 
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